Dual Order Parameter for the Nodal Liquid 
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The guiding conception of vortex-condensation-driven Mott insulating behavior is central to the 
theory of the nodal liquid. We amplify our earlier description of this idea and show how vortex 
condensation in 2D electronic systems is a natural extension of ID Mott insulating and 2D bosonic 
Mott insulating behavior. For vortices in an underlying superconducting pair field, there is an 
important distinction between the condensation of flux hc/2e and flux hc/e vortices. The former 
case leads to spin-charge confinement, exemplified by the band insulator and the charge-density- 
wave. In the latter case, spin and charge are liberated leading directly to a 2D Mott insulator 
exhibiting spin-charge separation. Possible upshots include not only the nodal liquid, but also a 
novel undoped antiferromagnetic insulating phase with gapped excitations exhibiting spin-charge 
separation. 
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I. INTRODUCTION 

The present paper is rooted in the conviction that the 
basic property of an insulator is that it insulate. Mag- 
netic order is a secondary effect which, though it may 
be one of the appurtenances of insulating behavior, is 
not synonymous with it. This premise underlies our re- 
cent discussion of the phase diagram of the underdoped 
cuprate superconductors, where we relied ppon the no- 
tion of an insulator as a vortex condensate^. We would 
like, here, to elucidate and expand upon this paradigm. 
As a by-product of this approach, we find a precise dis- 
tinction between spin-charge separated and confined in- 
sulators. 

Since correlated insulators often exhibit magnetism, 
they are typically described by their magnetic order pa- 
rameters. Furthermore, in commensurate, weak coupling 
models, the development of magnetic order can be the 
mechanism by which a system becomes insulating. This 
has led to the conflation of magnetic order with insulating 
behavior. However, this state of affairs is both unsatisfy- 
ing and incomplete since magnetic order can persist even 
when the system becomes conducting and, conversely, a 
system can be insulating even in the absence of mag- 
netic order, as exemplified by the nodal liquid. Hence, 
it would behoove us to seek an order parameter which is 
directly related to the electrical properties of an insula- 
tor. The difficulty in such a program is that an insulator 
seems more disordered than a conductor since most cor- 
relation functions in the charge sector are short-ranged. 
We interpret this as suggesting a 'dual' approach based 
on a 'disorder' parameter. Distilling the key elements of 
our nodal liquid construction, we propose that the ap- 



pearance of a non-zero expectation value for a 'disorder' 
parameter signals insulating behavior. 

The relation between the conduction properties of a 
system and spontaneous symmetry breaking in a dual 
order parameter is perhaps most transparent in the field- 
theoretical formulation of duality.cl In this version of the 
transformation, the dual theory is constructed to im- 
plement local charge continuity, which is a dynamical 
consequence of U(l) symmetry in the original Hamilto- 
nian, as a rigid constraint (dynamics in the dual theory 
implies conservation of vorticity). This construction is 
quite familiaii from the one-dimensional theory of Lut- 
tinger liquidsa. There, the two-current can be written in 
terms of a phase field (j>{x) as 
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(1.1) 



In the alternative, but equivalent, dual description, the 
current can be written in terms of a dual field, 6: 
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In a one-dimensional Mott insulator, the field acquires 
a mass, implying a gap in the spectrum of charge excita- 
tions, and hence, experimentally, in the optical conduc- 
tivity. Note that massive dynamics for the phase field <f> 
is inconsistent since it would violate charge conservation 
diji = 0. The one-dimensional model has thus commu- 
nicated an important lesson: insulating behavior occurs 
when a gap is acquired by the dual field representing the 
current operator. 

This description of Mott insulators generalizes readily 
to two-dimensional bosons. Let us write the boson anni- 
hilation operator, ip, in terms of its amplitude and phase: 
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-0 = ^/pe lip . In the superfluid state, we fix p and describe 
the system by its phase degree of freedom, <p(x). 



£ = \ (d^f ■ 



(1.3) 



We can model the destruction of superfluidity by intro- 
ducing a vortex field, <I>; effectively, we are reducing the 
amplitude degree of freedom to a vortex field which keeps 
track of the points at which-^t vanishes. We can represent 
the current, = d^ip, asB 



3n 



(1.4) 



1.2 to two- 



which is the natural extension of Eq 
dimensions. This parameterization is highly redundant 
as a result of its invariance under the gauge symmetry: 



a M -> a M + df+x, 



(1.5) 



for an arbitrary function, \. This gauge symmetry is 
enormously larger than the analogous global invariance 
in the one-dimensional Luttingcr liquid, 8^8 + const. 
The vortices see the gauge field, a M , according to the 
Magnus force law, so the dual Lagrangian takes the form 

C D = \ {d^f + - iaj 4>| 2 + V(9). (1.6) 

In the superfluid state, vortices are gapped and the 
bosons condense, while in the Mott insulating state, the 
vortices condense and the bosons become gapped (at the 
chemical potential). When the vortices condense, be- 
comes massive by the Anderson-Higgs mechanism. As a 
consequence, the system is insulating. 

Let us abstract away the archetypal featur es of this sys- 
tem. We can introduce the representation ( |1.4| ) for any 
conserved current in two-dimensions, so we can certainly 
use it for the electrical current in the fcrmionic system of 
our choice. In order to carry the rest of this scheme over 
to electronic systems, we must find a way for the system 
to conspire to make massive. The only way that this 
can hap pen which is consistent with the gauge symme- 
try (1.5), is via the Anderson-Higgs mechanism.u Since 
the Anderson-Higgs phenomenon can only take place if 
there is a condensate which is c oup led to a M according 



to the Magnus force law (Eq. F6), we are led to the 
following question: how do you define a vortex field in a 
fermionic system? One possibility is to implement statis- 
tical transmutation to represent the fermions as. bosons 
coupled to an auxiliary Chern-Simons gauge fielcB. Then, 
we can define vortices in the bosonic field. This approach 
is probably suitable for describing a conventional antifer- 
romagnet, as discussed very briefly in Section IV. But in 
this paper we pursue a different tack - using Cooper pairs 
as the bosons. This is quite promising for the cuprates 
because it is tailor-made for insulators which contain the 
germ of superconductivity. 

A new question rears its head when we consider an in- 
sulator which descends in this way from a superconduc- 
tor: do the finite-energy excitations inherit their quan- 
tum numbers from the superconducting state, or do they 



simply have the electron quantum numbers? In particu- 
lar, one can ask what is the energy of an isolated neutral 
spin-1/2 excitation. If this diverges with system size, 
then spin and charge are confined. If, on the other hand, 
it is finite, as it is in a superconductor, then the insu- 
lator exhibits spin-charge separation. A band insulator 
is, of course, a state of the former variety. As we show 
in the next section, it can be understood (rather differ- 
ently than in elementary textbooks) as resulting from 
the condensation of flux hc/2e vortices in a state with 
s-wave pairing. Spin and charge are confined as a result 
of the Aharonov-Bohm phase which a spinful excitation 
acquires as it orbits an hc/2e vortex; an isolated spin-1/2 
has logarithmically divergent energy in 2D. This spin- 
charge confinement physics is also present in a charge- 
density-wave (CDW) insulator, which occurs for exam- 
ple in the negative u (extended) Hubbard model at half- 
filling. There, however, the hc/2e vortex condensation 
leads to translational symmetry-breaking. In both in- 
stances the resulting insulating phase can be viewed as 
a "crystal" of charge 2e "Cooper" pairs. For 2D elec- 
tronic systems at or near half-filling with strong on-site 
repulsion, however, such CDW order is physically unrea- 
sonable, and for the cuprate materials can be discarded 
on phenomenological grounds. For this reason we are led 
to consider the possibility of the condensation of double 
strength hc/e vortices in a superconducting pair fielcB. 

The Mott insulator which arises at half-filling upon 
condensation of hc/e vortices has a number of appealing 
and remarkable properties. The insulating phase is trans- 
lationally invariant, in contrast to the CDW. Moreover, 
since there are no Aharonov-Bohm phase factors when 
a spin 1/2 is transported around an hc/e vortex, tho-.rc.- 
sulting Mott insulator exhibits spin- charge separationLTLl. 
For a d-wave superconductor appropriate to the cuprates, 
condensation of hc/e vortices leads directly to the nodal 
liquid. The nodal liquid indeed posseses gapless spin 1/2 
but charge neutral fermions - the nodons - which descend 
directly from the d-wave quasiparticles. As we shall see, 
there are also massive charge e spinless boson excitations 
in the nodal liquid phase. These "holons" are solitonic 
topological excitations in the underlying hc/e condensate, 
a dual analog of Abrikosov flux tubes. The excitations 
in the nodal liquid have the same quantum numbers as 
in the spin-charge separated gauge theoriest2l, but are 
weakly interacting, rather than strongly coupled by a 
gauge field. 

A peculiar feature of the nodal liquid is that spin- 
charge separation survives the ordering of the nodal spins 
into a phase with long-ranged antiferromagnetic order. 
This novel phase - denoted AF* - which has gapped 
nodons is distinct from the conventional Neel antifer- 
romagnet, AF, which does not have neutral, spin 1/2 
excitations even at high energiesEil. These two phases 
are physically very different, as may be seen from simple 
gedanken experiments which make the point that charge 
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e can be physically separated from spin 1/2 with finite 
energy cost in the AF* but not the AF phase. However, 
in two-dimensions nodon-holon bound states form in the 
AF* phase, so the spin-charge separation is not so easily 
found in the electron spectral function. 

Our ultimate goal is to describe a novel spin-charge 
separated state, the nodal liquid, and an ordered state 
which can result from it, AF*. Along the way, how- 
ever, we will re-examine a number of seemingly quotid- 
ian states such as the band insulator, the charge-density- 
wave, and the antiferromagnet. Our broader framework 
will enable us to understand the physics of doping these 
insulators from the point of view of cpating topological 
excitations in the 'disorder' parameterc3. Such a point of 
view naturally leads to a discussion of the possibility of 
spin-charge separation in these states. In sections II and 

IIIA, we will illustrate the physics of flux hc/2e condensa- 
tion in systems with attractive electron-electron interac- 
tions, where we expect the insulating states to be related 
to s-wave superconductivity. In the resulting states, the 
band insulator and the CDW, spin and charge are con- 
fined, as we discuss at length in section II. In section 

IIIB, we then consider flux hc/e vortex condensates in 
a d-wave superconductor, filling in a gap in our earlier 
paper. We introduce an effective lattice model in section 
IIIC which incorporates this physics. In section IV, we 
discuss the spin-charge separated antiferromagnet, AF*, 
and compare it to the conventional antiferromagnet, AF. 
Possible experimental signatures are analyzed. We con- 
clude with some summary remarks in section V, relegat- 
ing some supporting technical details to the appendix. 

II. THE BAND INSULATOR AND 
SPIN-CHARGE CONFINEMENT 

In the absence of electron interactions, a band insu- 
lator with two electrons per unit cell corresponds to a 
filled valence band of non-interacting levels. Provided 
the interactions are small compared to the energy gap to 
the conduction band this should provide a good descrip- 
tion of the phase. But even with stronger interactions 
a band insulator can be adiabatically deformed (without 
gap closure) back to the non-interacting state. To obtain 
an order parameter for the band insulator, we will at- 
tempt to describe this phase as a "quantum disordered" 
s-wave superconductor. 

To this end, consider spinful electrons moving in the 
two-dimensional continuum. In the presence of a local 
attractive interaction the Fermi surface is unstable, and 
presumed to form a spin-singlet s-wave superconducting 
phase, with energy gap A. Now imagine introducing a 
periodic potential with magnitude V and a period cor- 
responding to two electrons per unit cell. For V much 
smaller than A the superconducting phase should be 
stable, but with increasing V one expects an eventual 
quantum transition into a band insulator. This transi- 
tion can be described as a quantum vortex unbinding 



transition, analogous to the thermally driven Kosterlitz- 
Thoulcss transition. At zero temperature the vortices can 
condense, giving one an order parameter for the band in- 
sulator. 

Of particular interest is the nature of the gapped ex- 
citations in the two phases. As we shall discuss, in 
a superconductor Jl is possible to define charge neu- 
tral quasiparticlesE3, which carry spin 1/2. In an s- 
wave superconductor these "spinon" excitations are fully 
gappecH, but can be gapless (called nodons) in a d-wave 
superconductor which we return to in Section IV below. 
Since the Cooper pairs carry no spin in a singlet super- 
conductor, the Cooper pairs and spinons provide a nat- 
ural spin-charge separated description of the supercon- 
ducting phase. On the other hand, the excitations in a 
band insulator are electrons which of course carry both 
spin and charge - the band insulator does not exhibit 
spin-charge separation. We would like to try and un- 
derstand the mechanism whereby the separated spin and 
charge excitations in the superconductor become "con- 
fined" upon entering the band insulator. 

To address these issues it is convenient to consider a 
low energy effective theory for an s-wave superconductor 
in which both Cooper pairs and the gapped quasiparticle 
states near the Fermi surface are retained. The appro- 
priate Lagrangian takes the form: C = C c + L v + Ci, n t 
with 

C c = ci(x, T)[id t - V 2 - n]c(x, r), (2.1) 
£<p = y(M 2 , (2-2) 

C mt = | A|e^c T (x)cj. (x) + h.c. (2.3) 

Here c a denotes an electron with spin a, and ip is the 
phase of the pair field, with magnitude |A|. Integrating 
over high energy electron states well away from the Fermi 
surface, will generate dynamics for the phase field. The 
appropriate form of £ v at low energy is essentially de- 
termined by symmetry. Here we have retained the lead- 
ing order terms in a gradient expansion, with kq = k 
the compressibility and — Kj = v^kq a superfluid stiff- 
ness. Henceforth we will set the velocity v c = 1. In gen- 
eral a Berry's phase termel of the form £eerry = n-odtf 
is allowed (see section III below) but with one Cooper 
pair (two electrons) per unit cell of the periodic poten- 
tial no = 1 and the Berry's phase term can be dropped 
since exp(i J dtC^ c „ y ) = 1. 

Notice that the phase fluctuations are strongly coupled 
to the electron operators through <Cj n j. To decouple this 
interaction and to exhibit the spin-charge separation it 
is convenient to consider the change of variables: 

f a (x) =e l */\ a {x). (2.4) 

The Lagrangian becomes, 
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£ = C f + C v + ^J^ : (2.5) 

with 

£/ - &[idt ~ V 2 - M ]/ Q + IAI/tCi)/^) + ^.c, (2.6) 
and J M is a quasiparticle 3-current operator: 

Jo = fifa, J = ifiVf a + h.c. (2.7) 

This c urre nt is not conserved as a result of the anomalous 
term ( [2.3| ), but the spin currents, 

Jo = faUffi fl = tfKpVfp + h.c. (2.8) 

are conserved. Here we have assumed that the phase field 
is slowly varying and have dropped terms involving two 
spatial gradients of (p. 

The Lagrangian C f can be diagonalized as usual by a 
Bogoliubov transformation, and describes gapped quasi- 
particles. Since the / operators are electrically neutral 
but carry spin one-half, these excitations are "spinons" . 
The spinons are coupled to the phase fluctuations via a 
Doppler shift type term. In the superconducting phase 
these phase fluctuations are small, and will generate a 
weak interaction between the gapped spinon states. 

To quantum disorder the superconducting phase and 
arrive at a description of the band insulator, we will need 
to allow for vortices in the phase of the pair field. In two- 
dimensions vortices are simply whorls of current swirling 
around a core region. The circulation of such vortices is 
quantized, since upon encircling the core the phase <p can 
only change by integer multiples of 2tt. The "elemen- 
tary" vortices have a phase winding of ±2ir. Since the 
Cooper pairs have charge 2e, in the presence of an ap- 
plied magnetic field these vortices quantize flux in units 
of hc/2e. Inside the core of a vortex the magnitude of the 
complex order parameter |A| vanishes, but is essentially 
constant outside. Since the position of these "point-like" 
vortices can change with time, their dynamics requires 
a quantum mechanical description. Thus a collection of 
many vortices can be viewed as a many body system of 
"point-like" particles. Since positive (+1) and negative 
(— 1) circulation vortices can annihilate and disappear 
they behave as "relativistic" particles. There is a con- 
served vortex "charge" in this process, namely the total 
circulation, and an associated|-current. A duality trans- 
formation can be implementedenj in which the phase ip is 
replaced by a dual field, 9, which is the phase of a vortex 
complex field, $ ~ e 10 . In a Hamiltonian description, $ 
and $^ can be viewed as vortex quantum field operators 
- which destroy and create vortices. 

A crucial element in the duality transformation is the 
total electrical three-current, J^°* = K^d^ip+J^/2, which 
must be conserved even in the dual representation. This 
is achieved by expressing the current as a curl of a gauge 
field, a M B: 

4 ot = e^xduax, (2.9) 



which automatically implies the continuity equation 
d^J 1 " 1 = 0. This representation also introduces a gauge 
symmetry into the problem, — > + d^A. It is this 
gauge symmetry which is spontaneously broken in the 
band insulating state. 

In Ref. [l] the duality transformation was implemented 
in the presence of the Doppler-shift interaction between 
the Cooper pairs and the spinons, giving a dual La- 
grangian of the form: Cd = Cf + C v + C a with a vortex 
piece of the Ginzburg-Landau form, 

£vM = yl(^ - *27ra M )$| 2 - r|$| 2 - u\<P\\ (2.10) 

and 

C a = ^-{e 2 J ~b 2 ) + ^-J^ lluX dva x . (2.11) 
2«o 2ko 

Here e,- = (djao — doaj) and b — eijdiaj are dual "elec- 
tric" and "magnetic" fields. The dual magnetic field b is 
simply the total charge density (in units of the Cooper 
pair charge). The last term in C a is the only one cou- 
pling the spinons to the vortices. However, this dual La- 
grangian is not valid since it ignores a strong statistical 
gauge interaction between spinons and hc/2e vortices. To 
see this consider taking a spinon (/) around a closed loop 
which encircles an hc/2e vortex. Along this ci rcui t the 
phase ip winds by 2ir. Due to the 1/2 in Eqn. ( |2.4| ) this 
implies that the spinon / must change sign upon com- 
pleting this circuit. This can be formally implemented 
by introducing branch cuts emanating from each and ev- 
ery vortex, across which the fermion wavefunction must 
change sign. This represents a strong and long-ranged 
"statistical" interaction between the spinons and hc/2e 
vortices. 

The presence of this long-ranged interaction clearly 
invalidates the form of the dual Lagrangian. One is 
tempted to try and incorporate the branch cuts by intro- 
ducing two new gauge fields coupled together by a Chern- 
Simons interactions. The most natural way of doing this 
is to introduce a coupling J^a^ of the nodons to a gauge 
field which attaches half of a fictitious flux quantum 
to each vortex. However, such a coupling is not gauge 
invariant since the nodon current, J M is not conserved. 
To avoid this problem - but at the cost of breaking spin- 
rotational invariance - we couple a M to the z-component 
of the nodon spin current, which is conserved. Specifi- 
cally, let a M and a^ denote statistical gauge fields coupled 
to the spinons and vortices, respectively, with a modified 
dual Lagrangian: 

C D =£ f +£ a + C v {a^ + oj) + a M J* + C cs , (2.12) 

and a Chern-Simons interaction: 

C cs = 2a ll e l _ iv \d v a\. (2-13) 

The Chern-Simons term effectively attaches flux tubes 
with strength 1/2 to each of the spinons and vortices. 
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This follows from the equations of motion obtained from 
dCo/da^ = and dCo/da^ = which imply, respec- 
tively, 

^ = ^j;, (2.i4) 

and 

^ v xd v a\ = (2.15) 
Here, the vortex three-current is given by, 

j£ = Im[$*(^-ia M )$]. (2.16) 

Consider now trying to condense the hc/2e vortices. 
In the ground state the spinon (/-fermions) are gapped 
out with < J* >= 0, so one can presumably set a* = 0. 
Setting < $ >= $ corresponds to a spontaneous break- 
ing of the gauge symmetry, and leads to an effective 
Higgs Lagrangian: L = $ 2 ) K tl a 2 l /2. In terms of the 
dual Ginzburg-Landau theory this describes the "Meiss- 
ner state". But since the curl of a M corresponds to the 
total electrical current, this phase corresponds to an in- 
sulator - the band insulator - with a charge gap. 

If the dual Ginzburg-Landau theory is type II it will ex- 
hibit topological excitations corresponding to penetrat- 
ing quantized "dual" flux tubes. In the electronic insu- 
lator these correspond to gapped charge ±2e spin-zero 
states, which are two-electron bound states. 

But now consider an excited state in the insulator 
which carries spin 1/2. This can be created by adding a 
spinon at the origin by acting with /t (x = 0) . The pres- 
ence of a spinon induces a statistical gauge field from the 
Chern-Simons term: 

h s (x) = eijdial = - J = ^6 2 {x). (2.17) 

Since h s {x) corresponds to an applied "magnetic field" 
in the dual Ginzburg-Landau theory, adding a spinon is 
equivalent to the insertion of a solenoid carrying one-half 
of a (dual) flux quantum. Being in the Meissner state, the 
dual Ginzburg-Landau theory will tend to screen out this 
applied magnetic field by generating currents that induce 
an opposing internal field, b(x). This follows readily from 
the energy in the Meissner state which takes the form: 

E= ^1 J d 2 x\V9-2Tra-2na s \ 2 : (2.18) 

where we have put $ = $oe l9 . With both 9 = and 
dju = the energy in the presence of the solenoid will di- 
verge logarithmically with system size. Apparently, the 
energy of an isolated spinon diverges in the thermody- 
namic limit. But in the presence of an induced internal 
magnetic field the energy will be finite provided the inte- 
grated flux is precisely one-half of the dual flux quantum: 

/**(,) ~i. (2.19, 



In physical terms this corresponds to an induced cloud 
of electric charge with magnitude e. Evidently, in the in- 
sulating phase an isolated S z = 1/2 excitation will bind 
charge e to form a spin-up electron. The resulting exci- 
tation has finite energy Similarly, an isolated S z = —1/2 
will bind charge — e to form a spin-down hole with the 
same energy as the spin-up electron. This mechanism for 
confinement of spin and charge is reminiscent of the con- 
finement of charge and flux which occurs in the quantum 
Hall effect. In the bosonic Chern-Simons formulation of 
the quantum Hall effect, the confinement of charge and 
flux also arises via a Higgs mechanism - in this case when 
the composite Boson condenses. 

For the above Ginzburg-Landau/Chern-Simons theory 
there is another finite energy configuration in the pres- 
ence of a single S z = 1/2 - a 2tt winding in the phase 
9 of the vortex field together with an internal field of 
one-half quantum which aligns with the "applied flux" 
(J cPxb(x) = +1/2). This creates a finite energy excita- 
tion with spin one- half and charge — e, corresponding to 
a conventional spin-up hole. We can create a spin-down 
electron of equal energy in a similar manner. On physical 
grounds the energy of the S z = ±1/2 states should clearly 
be degenerate. Unfortunately, for the above Ginzburg- 
Landau/Chern-Simons theory while both energies are fi- 
nite, they will in general be different, due to differences 
in the core energies (eg. near the S z = —1/2 state with 
9 winding, it is necessary to suppress the magnitude of 
$ to zero). This signals a clear deficiency in the Chern- 
Simons formulation of the "statistical" interaction be- 
tween spinons and hc/2e vortices. Since the vortices and 
spinons can sense the sign of the statistical flux (that 
is ±1/2 flux quanta are not identical) the Chern-Simons 
fields do not give a faithful representation of the branch 
cuts. The necessary evil of breaking spin-rotational sym- 
metry is a consequence of this asymmetry in the Chern- 
Simons formulation. 

Currently, we do not have a convenient formulation 
of interacting hc/2e vortices and spinons which correctly 
respects this symmetry. Such a formulation would be 
particularly desirable for the case of a d-wave supercon- 
ductor where the quasiparticles are gapless. Neverthe- 
less, we believe that the Chern-Simons formulation does 
help elucidate the correct mechanism behind confinement 
of spin and charge upon condensation of hc/2e vortices. 

The preceding considerations dovetail naturally with 
the following approach to understanding spin-charge sep- 
aration. Consider the following gedanfeep-jexperiment, for 
which we are indebted to B. I. Halperinlii] , which probes 
the existence of spin-charge separation. Consider a to- 
tally gapped system which exhibits spin-charge separa- 
tion, meaning that it has weakly-coupled neutral spin 
1/2 excitations - to which we'll liberally apply the term 
nodons (even though the rest of our discussion applies 
also to systems which have nothing to do with d-wave 
superconductivity), following Ref. |l|- and charge e spin- 
less excitations - which we'll call holons, following Ref. M. 
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Imagine imposing two spatially localized perturbations. 
These take the form of an interaction Hamiltonian 



A 

+A 



(Q(x) 



(s z (x)y 



(Q(x')) 2 + (^(x')-l/2) 2 



(2.20) 



Here Q(x) = £ y + y) exp(-|y| 2 /^), and S 2 (x) = 
^ y s z (x + y) exp(— \y\ 2 /S, 2 ) are the total charge and z- 
component of the spin within a smoothed region of linear 
size £ around the point x. The perturbation favors local- 
izing a charge e without spin near x and a spin 1/2 with- 
out charge near x'. Now imagine taking |x — x'| 3> £ 3> a 
(the lattice spacing), so that the points are well sepa- 
rated. For small A, the ground state of the system will 
be unchanged, since there is a gap to all excitations. In- 
creasing A will ultimately induce a change in the ground 
state to take advantage of these perturbations. Provided 
£ is taken larger than the size of the nodon and holon, 
these excitations can come into the system to lower its 
energy and will be localized in the wells. This change in 
the ground state will occur at finite A, since the energy 
gap to the nodon and holon is finite, and indeed the crit- 
ical A will saturate as |x — x'| — » oo. One can interpret 
this critical A c as the minimum energy needed to produce 
an unbound nodon-holon pair. 

Now imagine repeating the same experiment on a band 
insulator or any other state which does not exhibit spin- 
charge separation. In this case, the nodon and holon are 
not available to "fill" the local perturbations. Instead the 
system must create a non-local superposition of elemen- 
tary excitations, i.e. develop a polarization, to localize 
the desired charge e or spin 1/2. Consider for example 
the region around the point x, in which a charge e should 
localize. A simple and generic model for a band insula- 
tor is a collection of deep potential wells, each containing 
two electrons. For the low-energy states, the wells may 
be approximated as quadratic and, in the ground state, 
all electrons are in the lowest harmonic oscillator level. 
The lowest excited states which do not induce local spin 
textures are constructed simply by moving both electrons 
in one of the wells from the ground state to one of the 
first excited states with energy gap loq. A convenient 
basis for these states is the set \fi) = aj„<4j0), where 

at = y/ mojo/2(x ll — ip^/muj) is the raising operator 
along the i axis in space. Now consider the superposi- 
tion of the ground and excited states, 



| U ) = (i + \uf/e 2 r 1/2 (lo) + 



■Im) 



(2.21) 



where t = \j2jmujn is the characteristic spatial width of 
the harmonic oscillator levels. For small \u\, the state |u) 
represents a small displacement of the electrons, i.e. 



uxu 



u. 



(2.22) 



Thus it possesses a dipole moment d = — 2eu, and a local 
charge may be built up by polarizing the collection of all 



the electrons near the point x (i.e. forming a many body 
state which is a direct product of single-particle states 
with differing displacements on each site). We consider 
a slowly- varying displacement u(x) for which a contin- 
uum description is adequate (although this is not a nec- 
essary restriction). Far away from the point x, the in- 
duced charge density p = 2eV • u, and hence the charge 



in a given region R is J R d d xp(x) 
radially symmetric configuration 



u(x") 



1 



2S d x - x 



2e J dR u ■ dii. The 



(2.23) 



thus carries a total charge e. Provided all the electrons 
are involved in the texture, there is no net spin polariza- 
tion. This polarized state is not an eigenstate of the un- 
perturbed Hamiltonian, but does couple favorably to the 
first term in Eq. 2.20. We can however determine the ex- 



pectation value of the (unperturbed) energy in this state. 
The result is essentially classical: E(u) = J d d xuJo(u/£) 2 . 
In two dimensions, this integral is logarithmically diver- 
gent, E(u) ~ InL. The isolated long-range polariza- 
tion thus cost an infinite energy. In the thought experi- 
ment, this divergence will be cut off by the finite distance 
between x and x', since we may localize an oppositely 
charged texture around the point x' in combination with 
an added electron with S z = 1/2, thus satisfying both 
perturbations and rendering the energy finite. However, 
the critical A c will grow logarithmically as |x — x'| — > oo, 
and hence it becomes impossible to create the isolated 
holon and nodon in the thermodynamic limit. 

This argument is appealing in that it agrees with the 
earlier Chern-Simons calculations which suggested loga- 
rithmic confinement of holons and spinons. As a means 
of distinguishing spin-charge separated from spin-charge 
confined phases, however, it is somewhat delicate. In 
particular, it fails for d > 2, where the polarization en- 
ergy to create the charge e texture becomes finite. It is 
also somewhat unsatisfying because the texture is not an 
eigenstate of the unperturbed Hamiltonian. Fortunately, 
the duality formalism allows the two phases to be distin- 
guished instead by the dual order parameter When 
<j> condenses, the statistical gauge interactions between 
nodons and the condensing vortices leads to spin-charge 
confinement in d — 2. In this sense $ is an order pa- 
rameter for confinement, a rather unique feature of the 
present theory. 



III. MOTT INSULATORS, HC/E CONDENSATES, 
AND THE NODAL LIQUID 



We now turn to the more interesting case of insulating 
phases with one electron per unit cell. In such Mott insu- 
lators electron interactions are necessary to destroy the 
metallic state, in contrast to the band insulator which has 
a smooth non-interacting limit. As in Section II we will 



G 



obtain a description of the insulating state by quantum 
disordering an appropriate superconducting phase. A 
brief discussion of quantum-disordered s-wave supercon- 
ductivity at this density illustrates the need to consider 
"double strength" vortex condensates in the cuprates. 
Such hc/e condensates are nodal Liquid insulators, and 
the subject of the remainder of this section. 



A. Quantum-disordered s-wave 

We begin by considering a system of spinful interact- 
ing electrons moving in the 2d continuum, which pair to 
form a spin singlet s-wave superconducting phase, and 
then "turn on" a periodic potential which for simplicity 
has square symmetry. Here, however, we choose the pe- 
riod to correspond to one electron per unit cell. As in 
Section II, with increasing potential strength the super- 
conducting phase can be destroyed by the unbinding and 
condensation of hc/2e superconducting vortices. But in 
this case there is only one-half of a Cooper pair per unit 
cell, so the resulting insulating phase will be dramatically 
different from the band insulator. In particular, one ex- 
pects the formation of a crystalline state of Cooper pairs 
which exhibits charge-density-wave ordering and sponta- 
neously breaks (discrete) translational symmetry. 

In order to proceed expeditiously to our main interest, 
the nodal liquid, we only mention a few salient points 
here. The difference between the CDW and the band in- 
sulator in the present approach is that the Berry's phase 
tcrmD, 



C-Berry = n Q dt^p, 



(3.1) 



cannot be dropped. To appreciate the physics of this 
term, we must return to the latt ice, where the lattice 
Hamiltonian corresponding to Eq. 2.2, together with the 
Berry's phase term (Eq. 3.1), takes the form: 



H v = — *2 ^ COs((^j - ipj) + U 2 *^2( n i - 



no) 



(3.2) 



(id) 



Here, n, is the Cooper pair number operator which is 
canonically conjugate to (p. Because, at half-filling, there 
is, on average, half a Cooper pair per site, one has 
n = 1/2. „ 

Implementing duality as beforeta (but now on the lat- 
tice) gives a d ual E uclidean acti on which is the lattice 
analogue of Eq. 2. 10| and Eq. 2.11, the only difference be- 



ing that the Maxwell term, ^-(ej — b 2 ), term is replaced 
by: 



S a = 



in 



(3.3) 



The new feature, as compared to the last section, is the 
"off-set" charge no which results from the Berry phase 
term ([O). It corresponds to an applied "magnetic" 



field for the lattice Ginzburg-Landau theory. The in- 
sulating CDW phase of the Cooper pairs corresponds to 
an Abrikosov flux lattice in this dual representation. As 
in the case of the band insulator, the vortex condensate 
leads to a charge gap and insulating behavior. 

Generally, one expects that a condensation of hc/2e 
vortices will lead to charge ordering with charge 2e per 
unit cell. This follows from the underlying duality: hc/2e 
vortices pick up a 2ir phase change upon encircling charge 
2e Cooper pairs - the same phase accumulated when a 
Cooper pair encircles such a vortex. Thus, the libera- 
tion and condensation of hc/2e vortices leads to a charge 
quantization in units of 2e. For the model with attrac- 
tive interactions and one electron per unit cell considered 
above, the resulting state is the (charge 2e) Cooper pair 
crystal. 



B. Quantum-disordered d-wave 

We now turn to the interesting problem of the 
quantum-disordered d-wave superconductor. D-wave su- 
perconductivity is thought to arise from the combination 
of strong on-site Coulomb repulsion and some unspecified 
(and controversial) longer-range attraction on the scale 
of a few lattice spacings. Certainly strong local repulsion 
is a key ingredient of the cuprates. For such systems, 
the hc/2e vortex condensation described above - which 
implies considerable double occupancy (at least in a re- 
gion near half-filling) - is physically unreasonable. It 
must also be discarded on phcnomenological grounds, as 
the actual undoped antiferromagnet is a Mott insulator 
without charge ordering. 

Although ft,c/2e-flux vortex unbinding is untenable for 
this case, phase coherence must nevertheless still be de- 
stroyed to obtain an insulating state. Charge uniformity 
and phase disruption can both be achieved together by 
unbinding bound pairs of vortices with flux hc/e instead 
of isolated ones. We expect such a double-vortex conden- 
sate to appear in the dual description as a condensate 
with a doubled dual charge, and hence a halved dual flux 
quantum!! At half-filling then, the dual lattice Ginzburg- 
Landau theory has a full 2tt flux per plaquette and thus, 
as desired, exhibits no translational symmetry breaking. 

Having motivated double- vortex condensation in the d- 
wave case, we now proceed to discuss its implementation. 
The calculations are significantly different because of an 
important additional physical ingredient in the d-wave 
superconductor: gapless fermionic quasiparticles. These 
excitations arise owing to the vanishing of the amplitude 
of the pair wavefunction at its nodes in momentum space. 
The presence of low-energy fermionic excitations necessi- 
tates a careful reinvestigation of the duality transforma- 
tion and its implications. Much of the necessary calcu- 
lations and formalism was discussed in detail in Ref. [I], 
and is briefly recapitulated in the appendix. 

As for the s-wave case, the analysis of the interactions 
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of vortices with quasiparticles is based on the neutralizing 
change of variables in Eq. 2.4. The distinctive feature of 
the d-wave superconductor is that the neutral spin- 1/2 
particles are gapless, and can be described by a Dirac 
Hamiltonian. Because in this case the spinons near the 
nodes can contribute to low-energy physics, we attribute 
to them special significance and the name nodons, sig- 
nifying the low-energy spinons descended from the d- 
wave nodal quasiparticles. Being gapless, they can be 
described by continuum field theory and a 4-component 
Dirac spinor ip (the analog of / in Sec. II - see the ap- 
pendix for a precise definition). 

Having already argued that only /ic/e-flux vortices 
should be considered, we will focus primarily on this sim- 
pler case. It is however, appropriate at this point to re- 
flect briefly on the consequences of these strong gauge 
interactions should single-strength vortices become im- 
portant low-energy excitations. As above, /ic/2e-vortex 
unbinding gives rise to strong gauge interactions of the 
spinons, now nodons. The gauge-theoretical arguments 
given in Sec. II can be again carried through, and we 
expect confinement of the nodons. Unlike the s-wave 
case, however, because the nodons are gapless excita- 
tions, their presence or absence has definite consequences 
on the ground state correlations. For instance, in a 
pure d-wave superconductor, the gapless Dirac excita- 
tions lead to static (equal time) power-law spin correla- 
tions and a T-linear magnetic susceptibility. It is natural 
to expect that the removal of the nodons from the low- 
energy spectrum in the insulator will be accompanied by 
the condensation of some pairing operator (for example, 
spin-density-wave order is characterized by the order pa- 
rameter (ij^ t v a y ail^) 7^ 0; other possibilities are legion), 
in most cases accompanied by a gap for the reconfincd 
electrons.E-3 The formation of such a paired-nodon state 
is analogous to chiral symmetry breaking in QCD, and in 
that context as well is generally believed to accompany 
confinement (although the converse need not be true). 
One can imagine approaching such a state from the su- 
perconductor by continuously lowering the vortex "mass" 
to zero, at which point one has a theory in which gapless 
fermionic nodons interact with gapless bosonic vortices 
via strong gauge interactions. This putative critical point 
is a tempting starting point for future systematic studies 
of such instabilities. 

We now return to the problem of he/ e vortex conden- 
sation, which, although motivated on energetic grounds, 
has dramatic consequences for the elementary excita- 
tions. Examining again Eq. 2.4 when only double- 
strength vortices are present, (p is defined modulo 47T, 
and this transformation defines a single-valued neutral 
fermion. The nodons thus experience no statistical gauge 
interactions in this case. The results obtained in Ref. [I], 
which ignored gauge interactions, hence apply to the hc/e 
condensate, with the proviso that the fundamental vortic- 
ity must be doubled throughout the analysis. The salient 
result is that when gauge effects are absent, the nodons 
and vortices interact only via the two-fluid interaction 



Lagrangian 



£-int — du'pJn 



(3.4) 



where J M is the electrical current carried by the quasi- 



particles, and is bilinear in the if) fields. Eq. 3.4 can be 



understood as the Doppler shift of the nodon energies in 
a superflow given by d^ip. This is a much weaker coupling 
than the statistical gauge interactions in the single- vortex 
condensate, and controlled analytic calculations are pos- 
sible. Detailed predictions for this quantum-disordered 
state, the nodal Liquid (NL), can be derived by writ- 
ing a coarse-grained continuum theory for $2 and o^, 
as in Refs. |l|,]l5|. The key conclusions are: (1) gapless 
nodons survive into the NL state, carrying spin but nei- 
ther charge nor current at low frequencies; (2) the NL has 
gapped charged excitations, the lowest-lying of which are 
expected to be (spinless) charge ±e holons, which occur 
as vortices in the dual order parameter $2; (3) the half- 
filled NL has uniform charge density, and upon (hole) 
doping charge is introduced as a spin-less Wigner crystal 
with charge e per unit cell (but see Sec. IV for a discussion 
of how this may be modified when antiferromagnetism is 
present). 

It is important to emphasize that a connection has 
been made here between two apparently unrelated phe- 
nomena. By assuming double- vortex condensation, char- 
acterized by the dual order parameters $2 = (e 2l0j ) 7^ 
(double vortices are condensed), $ = {e ') — (single 
vortices are bound), we were led to the persistence of 
spin-charge deconfincment in the insulator. The single- 
vortex disorder parameter, $ (which also distinguishes 
translational symmetry breaking at half-filling) can thus 
be regarded as an order parameter for conGnement. 



C. Lattice model for the nodal liquid 

We conclude this section by describing a direct route 
to the NL at half-filling, by which most of its properties 
may be derived without the use of the duality mapping. 
To do so, consider the following lattice regularization, 
which forbids /ic/2e-flux vortices from the outset: 

(ij) 

+ \A\(~l) x ^e~ l{lp ^ )/2 4- l c] l +h.c, (3.5) 
'(pi _ tpf 
2 2 



—j cos ( 



/ 2n» + c\ a c ia - I 



(3.6) 



The cos(ifii/2 — (fij/2) term has been chosen to allow 
±4-77 phase slips but not ±2ir phase slips, hence "con- 
fining" single-strength vortices. Further, we have made 
the apparently arbitrary choice of dividing the super- 
conducting pair-field phase amongst neighboring sites. 
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While this may seem unnatural, provided the contin- 
uum d-wave quasiparticle Hamiltonian is an adequate 
low-energy description, any lattice regularization should 
reproduce identical low- wavelength behavior. Finally, we 
have included a "charging energy" term coupling to the 
total (Cooper pair + quasiparticle) charge. 

This model has particle-hole symmetry, and at zero 
chemical potential is thus automatically at half-filling. 
To determine the properties of the system, we begin by 
performing the lattice analog of Eq. 2.4: cJ- Q = e lipj ^ 2 fj a - 
Simultaneously, to avoid non-trivial commutation rela- 
tions between /, /t and n, we let iVj = 2rij -I- cjc,j and 
<f>i = ifi/2. The / fermion creates neutral, spin-1/2 
quanta. In these variables, the Hamiltonian becomes 



H — H v + Hf + H,,, 



with 



U 



H ^ = ^ -J cos (fc - 

(ij) i 
(ij) 

H^^-t (e«+<-Mflf ja +h 
(ij) 



)+E^,-l) 2 , (3.7) 



(3.8) 
(3.9) 



Note that the nodon-phase coupling has been transferred 
from the pair-field interaction to the kinetic term by the 
operator transformation. An insulating state is obtained 
in the limit U 3> J, t, where the charging energy dom- 
inates over both pair and single-particle hopping. This 
state can be studied perturbatively in t and J, expanding 
around the insulating state with iVj = 1 exactly on each 
site. At t = J = 0, however, the /-particle Hamiltonian 
is still highly degenerate. This degeneracy is broken at 
second-order in t and J, giving the effective Hamiltonian 
(obtained, e.g., by perturbatively integrating out Ni and 
<j>i in a path-integral formulation) £T e // = H e o + H e i with 



(ij) 



V 



£iAK-ir-^c! T c]. i+ h, 

(ij) 

2t2 Vftf ft f 

/ j J iaJ jot J j(3J 0' 



u 



(3.10) 
(3.11) 



(ij) 



The quadratic Hamiltonian H £ q is identical to the mean- 
field Hamiltonian for d-wave quasiparticlcs, with a renor- 
malized bandwidth StJ/U. It therefore describes two sets 
of spin-1/2 Dirac fermions at low energies. One thus re- 
covers in this way the NL phase obtained previously via 
continuum duality. The interaction H e \ can be rewrit- 
ten as a combination of antiferromagnetic exchange and 
contact repulsion of the /-particles. If both are weak (as 
in the large U limit), such four- fermion interactions are 
strongly irrelevant around the non-interacting NL fixed 
point, due to the linearly vanishing density of states of 



the Dirac fermions. A slightly refined analysis includ- 
ing a physical external gauge field allows one to cal- 
culate explicitly the conductivity and show that the /- 
fermions caery no current in the NL state, so it is indeed 
an insulator!! 



IV. ANTIFERROMAGNETISM 

A. Phenomenology 

As discussed in the previous section, /ic/e-flux vortex 
condensation atop the d-wave superconductor yields the 
NL, an insulator with charge quantization (charge e per 
unit cell) appropriate near the half-filled Mott insula- 
tor. Unlike the CDW obtained by /ic/2e-flux vortex con- 
densation or the conventional (fully gapped "short-range 
RVB" ) spin liquid state,! the NL also contains low-lying 
gapless spin degrees of freedom, the nodons, which con- 
tain the germ of true antiferromagnetic order. As de- 
scribed in Ref. |l|, the phenomenology described above 
can be easily extended to include Noel order. 

The description to this point has essentially neglected 
inter-nodon interactions. This approach is justified pro- 
vided such interactions are weak, as all such terms are 
perturbatively irrelevant (in the renormalization group 
sense) in the Dirac theory describing the NL.El However, 
perturbative irrelevance does not imply that strong in- 
teractions cannot drive quantum phase transitions and 
hence a qualitative change in behavior. Indeed, interact- 
ing Dirac fermion models are known to undergo chiral 
symmetry*-hreaking transitions as quartic couplings are 
increased.!!^ The nature of the transition incurred de- 
pends upon the precise nature of the interactions, and 
various circumstances can induce antiferromagnetism, 
spin-Pcicrls, charge density wave, and other types of or- 
dering from the NL Lagrangian. Indeed, in the lattice 

L) capa- 



model above we obtained an interaction (Eq. 3.11 



ble of driving a transition to an antiferromagnetic state 
if sufficiently large. 

Because of the uncertainties and pitfalls of attempt- 
ing a microscopic justification of such interactions, how- 
ever, we prefer to follow the strategy of Ref. [I] and take 
a more phenomenological approach. For simplicity let us 
focus on the case of half- filling with particle/hole sym- 
metry. Since the cuprate materials are clearly antifer- 
romagnetically ordered, we will assume the existence of 
a triplet collective mode with momentum (tt,tt). In the 
NL phase this magnon mode has a gap, and indeed we 
expect also a non-zero lifetime, as the triplet magnons 
can decay into pairs of nodons, so strictly speaking the 
magnons are not sharply defined elementary excitations 
in this case. Nevertheless, we may imagine tuning a pa- 
rameter (e.g. reducing frustrating spin-spin interactions 
in a lattice model) to reduce the magnon gap. Ultimately 
when it vanishes the collective mode becomes sharp and 
condenses to form an antiferromagnetically ordered state. 







In the antiferromagnetic state, the non-zero Neel vector 
coherently mixes nodons with opposite quasi-momentum 
and opposite spin, halving the magnetic Brillouin zone as 
is usual in spin-density-wave systems. As in those more 
conventional cases, this has the effect of opening up a 
gap in the nodon spectrum. 



E n (q) = ±^(v F q ± r + (v A q^ + (gN o r, (4.1) 

where vf and va are the Dirac "velocities" perpendicu- 
lar and parallel to the putative Fermi surface, q is the 
momentum measured from a node, Nq is the mean-field 
staggered magnetization, and g is a phenomenological 
coupling constant. Due to this mixing, the only gapless 
degrees of freedom in the antiferromagnet are the col- 
lective spin-wave modes guaranteed by Goldstone's the- 
orem. 

In this way we arrive at an effective low-energy field 
theory for an antiferromagnetic Mott insulator, which we 
will denote (for reasons which will become apparent) as 
an AF* phase. The antiferromagnetic spin order, feature- 
less incompressible charge configuration, and gapless spin 
waves are qualitatively identical to those we would obtain 
in more conventional antiferromagnetic models, e.g. the 
nested spin-density-wave in the weakly-interacting half- 
filled Hubbard model, or alternatively the t — J-like very 
large U limit of the same Hamiltonian. We stress, how- 
ever, that although the nodons have been lifted away 
from zero energy, they are not confined by the sponta- 
neous symmetry breaking in the AF* phase, i.e. spin- 
1/2 neutral particles still exist as well-defined elemen- 
tary excitations. For this reason, we believe that such a 
novel antiferromagnetic insulator (AF*) is topologically 
distinct from (i.e. cannot be adiabatically deformed into) 
a more conventional antiferromagnetic (AF) state. This 
conviction is bolstered by the existence of the dual or- 
der parameter <E>2, which we have argued characterizes 
the nodal liquid and the AF* phase. Since $2 creates 
hc/e vortices in a pair field, the AF* phase contains the 
germ of superconductivity. In contrast, construction of a 
dual order parameter for the conventional antiferromag- 
net probably requires the use of Chern-Simons (charge 
e) bosons, obtained directly by statistical transmutation 
from the electrons. For example, condensing elementary 
hc/e vortices in the spin-up boson to form a charge e crys- 
tal which lives on one sub-lattice, and similarly freezing 
the spin down particles onto the other sub-lattice, should 
suffice to describe the hidden order of the AF phase. Ev- 
idently, the dual order parameters in AF and AF* are 
very different. Alternatively, one can distinguish AF* 
and AF by testing for the presence or absence of spin- 
charge separation employing the argument in Section II. 
To make the argument precise in this case probably re- 
quires adding an easy-axis anisotropy, 

H ea = J2Jea[(s X ) 2 + (s V f]- (4-2) 
x 

which creates a gap in the magnon spectrum. But it is of 



more interest to address whether the AF and AF* phases 
can be distinguished experimentally? To address this, we 
turn to a discussion of the electron spectral function in 
these two phases. 



B. electron spectral function 

In standard many-body systems, the existence of well- 
defined excitations is ascertained by examination of 
the relevant spectral function.EZI Unfortunately, a direct 
probe of spin-charge separation via spectral functions is 
not possible, since there are no local operators which sep- 
arately create nodons and holons. On the other hand, the 
electron spectral function, A(k,uj), is accessible experi- 
mentally and has been intensively studied in the high- 
temperature superconductors with momentum resolution 
via angle-resolved photo-emission spectroscopyta~E3 and 
locally (i.e. in momentum-integrated form) via non- 
linear tunneling characteristics.Eillt seems natural to sug- 
gest that A(k, bj) might possibly give one a way to distin- 
guish the AF and AF* phases, since in the latter unbound 
nodon-holon pairs form a two-particle continuum, and in 
the former the electron is itself the elementary excitation. 



1. AF 

To see how this idea works out in practice, let us first 
consider in some detail the spectral function in the AF 
phase. A simple model which captures the qualitative 
physics of the spectral function is fluctuation-corrected 
spin-density-wave mean field theory. The (imaginary 
time) quasiparticle Lagrangian is 



L 



4 [d T +e k - l-i] c k + gN 



C \L+Q at 



k- 



(4.3) 



where we have assumed ordering at Q = (jr, 7r). For 
simplicity, let us assume the Fermi surface intersects 
(7r/2,7r/2) with some curvature (it is straightforward to 
generalize this to other geometries). Choosing new coor- 
dinates along the (1, 1) and (1, —1) axes, we then write 
k = (n/y/2, 0) + q, and £k — M ~ V F 1x + ?y/2m near this 
point. Performing a similar expansion near the opposite 
point on the Fermi surface, and defining continuum fields 
f/ ± (q) = c (±7r/y 2 0)+q , one finds 



C-Nqp = fy 



d T + ivFT z d x 
+gN ■ tfaT x fj, 



2m d y 



(4.4) 



where g is the spin-density- wave coupling constant. 
We suppress the additional quasiparticles located near 
(±7r/2, — 7r/2), since these are not coupled to the 17 fields 
by the ordering wavevector Q. In the antiferromagnetic 
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phase, < N >= Nq 7^ 0, and if fluctuations are ignored 
the electron states are gapped with energy dispersion, 



q 2 y /2m, 



(4.5) 



with A = gNo the mean field spin-density wave gap. 

Spatial and temporal fluctuations of the Neel field 
N can be d escribed by, e.g. a Landau theory such 
as Eq. A10, or by a nonlinear sigma model. In the 



AF phase, we require only the spin- wave expansion for 
small deviations, IL <C 1, from perfect alignment, i.e. 
N = iVo(ni,n 2 ,Vl - n 2 ), for small U. Since uniform 
rotations of N are equivalent by SU(2) invariance, it is 
convenient to perform the "gauge" 



rotation 



ry(x, r) = exp[iei J -cr i II ;; -(x ) r)]»7(x,r). 



(4.6) 



In the new variables, the quasiparticle Lagrangian be- 
comes LNqp 



C n + Ca- V , with 



I + Aa z T a 



-leijT)' 



dyU-idy 



(4.7) 
(4.8) 



governed by the quadratic Lagrangian 



Here / = d T + ivpr z d x — ^dy. Finally, the magnons are 



K 
~2~ 



|a r n| : 



<ivni 



(4.9) 



The spin- waves in Eq. 4.9 are gapless as required by Gold- 
stone's theorem. Neglecting the coupling to the spin- 
waves (a good approximation if K is very large, so that 
the IT fields fluctuate very little), the r\ particles are non- 
interacting quasiparticles with a gap A, and have a sharp 
spectral function 

A°(q, u) = TT^ImG^q, iw n -> w + iS) = S(uj - E(q)). 

(4.10) 



The spin- wave coupling, Eq. 4.8 , generates a self-energy 
in the rj Green's function. The relevant diagram is indi- 
cated in Fig. 1, and a straightforward if tedious evalua- 
tion shows that 



G 7) (q, uj) = [iuJ + v F T z q x + q 2 y /2m 
+Aa z T* + S(q, W )]- 1 , 

ImE(q,c i ;)~-4 T T(H d 0(M 5 

XV Da 



(4.11) 
(4.12) 



with 5u = uj — E(q). Eq. 4.12 holds provided \q\ < mv s . 
Since ImE <C 6u> for small 6u>, the decay rate is negligible 
at low energies and we expect a delta-function singularity 
to survive in the spectral function at Su> = (there will 
also be some small shifts in the energy spectrum itself 
given by the real part of S). 



Fig. 1: Self-energy diagram for the charged quasipar- 
ticles in the AF phase. The solid line indicates a dressed 
?7-fermion Green's function, while the dashed line is the 
spin-wave propagator. 



When calculating the electron spectral function, one 
needs to include the effects of the SU(2) rotation, Eq. 16 . 
Since the physical electrons are created by the r) fields, 
additional factors of the II operators appear in A(k, oj). 
The LI operators may be expanded out of the exponen- 
tial and treated perturbatively. Both this effect and the 
broadening due to the non-zero self-energy above lead to 
additional weight for Slu > 0, which is often referred to 
as "incoherent" spectral weight. The general expecta- 
tion for the electron spectral function in the AF phase is 
illustrated in Fig. 1. 

The physical meaning of these results is the follow- 
ing. The minimum energy excitation with charge e and 
spin-1/2 is the electron, which in the interacting system 
is "dressed" by magnon excitations that mix with the 
bare electron. Further, there are higher-energy states in- 
volving a dressed electron and unbound excited magnons 
which arc orthogonal to the interacting electron but not 
the bare one, and thus show up as continua for Slu > 
once interactions are present. The true elementary exci- 
tation does not decay however, basically because phase 
space (which leads to the ImE ~ (Suj) d law above) pre- 
vents it. Thus the expected electron spectral function has 
a resolution-limited dispersing peak at the single-particle 
gap near its minimum, above which lies continuous spec- 
tral weight. Well away from (tt/2,tt/2), phase space may 
(or may not) open up to allow decay even of the single- 
particle peak, depending upon details of the band struc- 
ture and interactions. 



2. AF* 

In the unconventional antiferromagnet, we expect the 
presence of unconfined nodons and holons to lead to a 
two-particle continuum in the electron spectral function. 
While this is indeed the case, lower-energy features in 
fact exist due to nodon-holon bound states. Such bound 
states are analogous to excitons in semiconductors, which 
provide sharp peaks in the optical conductivity despite 
the existence of an electron-hole continuum at higher en- 
ergies. 

Similar considerations apply here. In particular, if we 
consider the interaction of a charge e holon with a spin- 
1/2 nodon, it is quite natural to expect that they may 
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experience an attractive interaction leading to a bound 
state with both charge e and spin-1/2, i.e. an electron. 
To show the existence of such bound states we specialize 
again to the case of half-filling with particle-hole symme- 
try. It is convenient to perform a particle-hole (Bogoli- 
ubov) transformation on the nodon operators near one 
pair of nodes, 

= V'lat; i>ai = i>\ al , (4.13) 

so that the number operator of the transformed fermions 
is proportional to the z-component of spin: S z = 
(1/2) J cPxip^ip. In the presence of Neel order, N = Noz, 
the transformed nodon Hamiltonian density takes the 
simple form, 

U n = ftHj, (4.14) 

with a single particle Hamiltonian: 

ffn(pi) = v{r z pl + r xP \) + A„rV». (4.15) 

Here, the nodon momentum operator pi = — iV ri is con- 
jugate to the position ri, and for simplicity we have as- 
sumed only a single nodon velocity. This Hamiltonian de- 
scribes massive nodon states, with energy gap A„ = gNo. 
Since the holons are also gapped, the appropriate first 
quantized Hamiltonian for a single holon (with position 
r 2 and momentum p 2 ) is simply 

P 2 

H h (p 2 ) = A h + ^—. (4.16) 

The form of the interaction between the nodons and 
holons follows from the dual Lagrangian in Section 2. 
For simplicity we only retain the density-density interac- 
tion term, proportional to J^jdiaj, where eijdiCij is the 
holon density and the nodon density can be expressed 
in terms of the transformed fermions as: J = rjp t z a z : ip ■ 
The corresponding first quantized interaction Hamilto- 
nian is then, 

H int ( ri - r 2 ) = ua 2 T z a z S^(r 1 - r 2 ), (4.17) 

with interaction strength u, and a a short distance cutoff. 

Since the two-body Hamiltonian is independent of the 
"center of mass" coordinate R = (ri + r 2 )/2, the total 
momentum, P = pi + p 2 , is conserved. For simplicity 
we consider bound states with P = 0. The Hamiltonian 
for the relative coordinates, 

r = n-r 2 ; p = ( Pl +p 2 )/2, (4.18) 

then takes the simple form 

H rel = H n (p) + H h (p) + H mt (r). (4.19) 

To solve for bound states with energy E, we recast the 
Schrodinger equation H re i<j) = E<j> in the form: 

G _1 (q)0(q) = -ua 2 T z a z (f>(r = 0), (4.20) 



with matrix Greens function G _1 (q) = H n (q) + £f/,(q) — 
E. Here </>(q) denotes the Fourier transform of the 
4— component wave function 4>(r). Upon matrix inver- 
sion this can be rewritten as M</>(r = 0) = 0, with 

M = 1 + ua 2 J ^L G (q)rV, (4.21) 

so that the eigenvalue condition is the vanishing of the 
determinant: det(M.) = 0. Here we are implicitly assum- 
ing that the integration is cut off at high momentum by 
q c = l/a. 

An explicit expression for the bound state energy, Eb, 
can be readily obtained in the u — > limit by putting 
Eb = A„ + A/, — £b with small binding energy, eb- In this 
limit one need only retain the contribution to the above 
integral which is infrared divergent, which gives 

M = 1 + W(<t z t z -<7 x t x ), (4.22) 

with 

W = -(u/8ne )He b /e Q ). (4.23) 

Here we have defined an energy scale, 

e = {m h v 2 + A n )/(2m h A n a 2 ). (4.24) 

Since det(M.) = 1 — AW 2 , the eigenvalue condition re- 
duces to W = 1/2, which gives the final result for the 
bound state binding energy, 

e b = e exp(-47re /u). (4-25) 

Notice that the binding energy is exponentially small in 
the interaction strength u, reflecting the two-dimensional 
constant density of states for free massive holons and 
nodons. If one were to change the sign of the interac- 
tion there is still a bound state (from W = —1/2) with 
the same energy. In either case, the bound state has 
the quantum numbers of the electron with s z = 1/2 and 
charge ±e. A spin down bound state can also be readily 
found, corresponding to the binding of a holon to a single 
nodon "hole" in the filled Fermi sea. 

Between the threshold energy for generating the elec- 
tron (Ef,) and the energy of the unbound nodon-holon 
continuum, the electron spectral function should be gov- 
erned by qualitatively the same physics as in the AF case, 
except that the total spectral weight of the correspond- 
ing feature will be reduced by matrix element factors 
arising from, e.g. the possibly large spatial extent of the 
nodon-holon bound state wavefunction. Upon reaching 
the nodon-holon continuum, we expect a much enhanced 
spectral weight but no sharp feature at the continuum, as 
it is already lying in the continuum formed by the elec- 
tron plus spin-wave excitations, and the pair excitation 
can thereby easily decay. 

The A(k, uj) in the AF and AF* phases are thus not 
qualitatively different, and cannot strictly speaking be 
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used to distinguish the phases. Quantitatively, however, 
we expect the AF* spectral function to exhibit a very 
small "quasiparticle" peak, with minimal separation from 
a nodon-spinon continuum carrying most of the spectral 
weight. If we assume that the holon gap greatly exceeds 
the gap for the nodons, then both features are expected 
to disperse in approximately d-wave fashion, though the 
cusp for angles near ±45° should be rounded by the 
nodon gap (Eq. p~l| ). At non-zero temperatures, ther- 
mally excited particles will scatter the injected electron 
and lead to broadening of even the threshold peak. In the 
AF* case, where this feature is expected to lie close to 
the nodon-holon continuum and have small weight, such 
thermal broadening could well remove the quasiparticle 
peak completely at experimental temperatures. 

For a system at half-filling but without particle-hole 
symmetry, the Neel ordering wavevector is not commen- 
surate with the spacing between antipodal nodes. If this 
incommensurability is sufficiently large, it is possible for 
the nodons to remain gapless even in the presence of long- 
range antiferromagnetic order. In this unusual state, 
which we denote as AF/NL, gapless spin 1/2 nodons 
co-exist with the spin 1 magnons. Since the density of 
states for the gapless nodons vanishes linearly with en- 
ergy, a weak interaction with the massive holons is not 
expected to result in a holon-nodon bound state. Angle 
resolved photo-emission in the AF/NL phase will thus 
have a number of notable features. Specifically, since the 
electron will decay into the nodon-holon continuum, one 
does not expect any sharp features in the momentum 
resolved spectral function. The lowest energy spectral 
weight is expected at the nodes, with a threshold energy 
which disperses linearly away from the nodes as in the 
d-wave superconductor. The spectral weight should rise 
smoothly above threshold due to the nodon-holon con- 
tinuum - with no delta function peaks. This behavior 
is in fact reminiscent of that observed in the undoped 
Ca compound by Shen et. al.£3 and is in marked con- 
trast to the delta-function spectral features expected in 
a conventional antiferromagnetic. 



V. DISCUSSION 



conductor in any dimensionll! For the superconductor 
the mechanism is different: Pairing of electrons into sin- 
glets creates a gapless collective (second) sound mode 
that carries the charge. The sound mode can adjust al- 
most instantaneously to a quasiparticle, effectively neu- 
tralizing it, leaving only a bare spin-1/2. On the face of 
it this species of spin-charge separation appears consid- 
erably different from the other topological varieties. 



Type I 



Type II 





Fig. 2: Phenomcnological phase diagrams in the nodal 
liquid Ginzburg-Landau theory of Ref. [l| Here x is the 
hole doping and T is temperature. As discussed in the 
text, hc/e vortex condensation leads to the unconven- 
tional AF* antiferromagnetic Mott insulator. Note that, 
depending upon the magnitude of particle/hole asymme- 
try, nodons may remain gapless in an AF*/NL state at 
half-filling - space constraints prevent us from indicating 
this on the figure. Depending upon microscopic parame- 
ters, two principle phase diagrams occur opon doping. In 
the type I scenario, added charge segregates into locally 
superconducting regions, which coalesce at some critical 
doping x c . In the type II scenario, added charges order 
into a Wigner Crystal (WC) with charge e per period, 
presumably with some associated spin ordering. After 
a small amount of doping the antiferromagnetic order is 
suppressed and the nodons are liberated into a nodal Liq- 
uid (NL) coexisting with the WC. At x c this WC melts 
into the d-wave Superconductor (dSC) phase. See Ref. |l| 
for details. Some modifications are necessary if the effects 
of impurities are included, some of which are discussed 
in Ref. M 



A small number of examples of condensed matter sys- 
tems are generally agreed to exhibit exotic quantum num- 
bers, i.e. particles which seem to require "splitting" the 
electron. Both charge fractionalizationEI and spin-charge 
separation are generic in one dimensional. In the two- 
dimensional quantum Hall effect, fractionally charged, 
particles have been known to exist for some time,E2l 
and recentlyJiave been observed in dramatic shot-noise 
experiments.E-3 In both these examples, fractional charge 
is connected to topological excitations: solitons or do- 
main walls in one dimension and vortices in two dimen- 
sions. 

A third example, less widely appreciated, is a super- 



In this paper we have exploited a dual formulation 
to show that indeed isolated charges ("holons") derived 
from the superconductor can be understood as topologi- 
cal excitations in a vortex condensate. Further, we have 
described how spin-charge separation can occur in an 
insulating state which results from the quantum disor- 
dering of a superconductor. Of course, propinquity to 
the superconducting state doesn't guarantee the inheri- 
tance of spin-charge separation: it only occurs when flux 
hc/e vortices condense. The most interesting example 
of this phenomenon - from the point of view of high-T c 
phenomenology - is the nodal liquid, which we have dis- 
cussed from this standpoint. The condensation of hc/2e 



13 



vortices, on the other hand, leads to the confinement of 
spin and charge. The band insulator and the CDW, for 
example, can be understood in this way. 

One striking consequence of the distinction between 
spin-charge separated and confined systems is that there 
are two distinct antiferromagnetic states: one, AF, which 
is the ordinary antiferromagnet and another, AF*, which 
is spin-charge separated. The latter results from or- 
dering the nodons in a nodal liquid, so it has neutral, 
spin 1/2 excitations. However, the distinction between 
the AF and AF* phases is experimentally rather elu- 
sive. The natural place to look is the electron spectral 
function, which can be probed through angle-resolved 
photo-emission experiments. Under some conditions (see 
Sec. IVB), the unconventional antiferromagnet should ex- 
hibit only a nodon-holon continuum instead of a quasi- 
particle pole. Unfortunately, the existence of nodon- 
holon bound states makes the distinction between the 
AF and AF* phases rather subtle. On the other hand, 
according to the paradigm presented here, spin-charge 
separation is a consequence of hc/e vortex condensa- 
tion. Thus, spin-charge separation could be indirectly 
evidenced by the observation of hc/e vortices near the 
quantum critical point at which superconductivity is de- 
stroyed. 

The continuity of spin-charge separation that is em- 
bodied in the nodal liquid and its offspring AF* state 
makes possible a simple phenomenological description 
of the evolution from the insulator to the superconduc- 
tor, as espoused in Ref. [l]. For example, the simplest 
Ginzburg-Landau formulation predicts the phase dia- 
grams in Fig. 2. If, on the other hand, spin-charge sepa- 
ration is absent in the undoped insulator, there must be a 
confinement transition between x = and the supercon- 
ductor. While we have argued that such a transition is 
driven by hc/2e vortices, its nature and the phases which 
it connects are highly non-trivial. The relative simplicity 
and elegance of the nodal liquid scenario thus argues in 
favor of its relevance to the cuprates. Despite numerous 
and interesting differences among different compounds, 
the phase diagrams of high-temperature superconduc- 
tors enjoy a remarkable degree of universality. A num- 
ber of theoretical works have attempted to understand 
the commonalities and variations amongst the_topology 
of these phase diagrams phenomenologicallyEj At low 
temperatures, however, we believe classical phenomenol- 
ogy based only on conventional order parameters misses 
the important physics of gapless quasiparticles and spin- 
charge separation that are key in the vicinity of d-wave 
superconductor-insulator transitions. Any viable theory 
of the cuprates must at least address the issue of how spin 
and charge either remain separated or become confined 
on approaching the insulator. 

In addition to the search for a 'smoking gun' experi- 
ment for spin-charge separation, there are a number of 
other interesting questions raised by this work. How do 
we implement the interaction between nodons and hc/2e 
vortices in an 5?7(2)-invariant way? How do these formu- 



lations of spin-charge separation apply to 31? systems? 
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APPENDIX: DIRAC LAGRANGIAN 

Here we review the effective Lagrangian for low-energy 
d-wave quasiparticles, following the notation of Ref. [y. 
It is most directly written in terms of the appropriate 
Nambu-Gorkov-like spinor, ^i aa , with 



*ii Q (k) 
* j2Q (k) 



°Ki+k,a' 

■ y t 

lcr Q( 3 c _(K I +k)^' 



(Al) 
(A2) 



where Ki and K2 are the momenta of the d-wave nodes 
along the Fermi surface. We use index-free notation in 
which Pauli matrices /t, r, a act in the node, particle-hole, 
and spin (iaa) subspaces, respectively; furthermore, if a 
single index is given explicitly, it is always the node in- 
dex. In a particle-hole symmetric model at half-filling, 
Kx/2 = (±7r/2,7r/2) in the usual (a, 6) crystalline co- 
ordinate system (i.e. axes along the Cu-0 bonds). The 
separation in Eq. |A2| is well-defined provided the momen- 
tum is restricted to points near the nodes, i.e. |k| < A, 
where A is a cut-off. 

As in the s-wave case, we must allow for space-time de- 
pendence of the superconducting phase (p. For a d x 2_ y 2 
superconductor, one has 

v(cii{t)c n {t))' = A d ( Xi -Xj)exp[i(p(x,t)], (A3) 

where the prime on the angular brackets indicates an av- 
erage (path integral) over high-energy electronic states 
away from the nodes, and x = (x; + Xj)/2. The ampli- 
tude function Ad(x) is the Fourier transform of the usual 
momentum-space gap function, Ak ~ /(|fc|)[cos 2 k a — 
cos 2 fcf,], and decays on the scale of £. It is usually 
more convenient for us to work in rotated coordinates 
x = (x a + Xb)/V2, y — (xb — x a )/y/2. The appropriate 
effective quasiparticle Lagrangian density was derived in 
Ref. § 

£ * = Yj *i ( idt + iv F T " d x + iv & T s e ls ^/ 2 d y e ls ^/ 2 ^ 9 t 

s=± 

+ (l^2;x^y). (A4) 
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Eq. A4 is derived on the assumption that the phase ip 
is slowly varying on the scale of the coherence length, 
i.e. \£d x <p\,\£d y <p\,\£/vd T tp\ <C 2n. However, we expect 
on grounds of universality that Eq. A4 and its conse- 
quences provide a correct low-energy description of the 
d-wave SG-and its quantum-disordered descendents more 
generallyEO 

The analysis of the interactions of vortices with quasi- 
particles is based on the important change of variables 



tp = exp(-i(^r z /2)*. 



(A5) 



Inserting Eq. |A5j into Eq. |A4j, one finds = + Ci nt , 
with 



Cj, = ip\[idt + v F r z id x + v A T x id y ]ipi 
+ (1 <-> 2, x «-> y). 



(A6) 



The nodon field ip i nter acts with the phase of the order- 
parameter as in Eq. 3.4. Here the electrical 3-current J M 
is given by 



Jo = 



(A7) 



V N {\N\) = r N \N\ 2 + u N \N\ 4 . The parameter tat controls 
the presence or absence of AF order. In mean-field the- 
ory, and neglecting for the moment the nodon coupling 
g, the ground state passes from long-range to short-range 
AF order as rjv is tuned from negative to positive. We 
include only the most relevant coupling of the Neel field 
to the nodons allowed by symmetry, 



C 



nodon 



= c„ 



gN ■ Sff, 



(A12) 



with N = (N). 

A compelling feature of the above description is the re- 
sulting low- lying spectrum in the antiferromagnet. The 
model can be readliy diagonalized with an appropriate 
Bogoliubo v tr ansformation, giving the energy eigenval- 
In all nodon sectors there is a non-zero 
The nodons having been lifted to 



ues,in Eq. 4.1. 
gap, equal to gNo. 
finite energy, the only remaining gapless excitations in 
the AF* phase are the spin waves (slow rotations of N) 
dictated by Goldstone's theorem. 



■ VF /t / 



(A8) 



Compared to the statis tica l gauge interaction with 



±/ic/2e-flux vortices, Eq. 3.4 represents a much weaker 
two-fluid interaction between the quasiparticle or nodon 
current J M and the supcrfluid current d^ip. A continuum 
duality transformation appropriate for such a coupling 
was described in detail in Rcf. [l], and on the lattice in 
Rcf. ISL Noting that the (Euclidian) nodon current J M 
couples to the superfluid current d^ip in a manner di- 
rectly generalizing the Berry's phase coupling inop, the 
dual (Euclidian) lattice action can be determined simply 



by replacing in^S 



(j,Q 



inoS 



fj,0 



S a — ► S a — ~y^ y (e^\A v c 



Nq5 



f i0 



(A9) 



The incorporation of antiferromagnetism was also de- 
scribed in Ref. [j]. A low-energy effective Lagrangian de- 
scribing the magnon mode and its coupling to the nodons 
is 

£ = ^K^N\ 2 - V N (\N\) + gN ■ (A10) 

where K = K, K\ = K 2 = —v 2 K, with v s the spin-wave 
velocity in the AF. Here 



s--i 



(All) 



is the spin operator at momentum 7?. Near any phase 
transitions, and for most phenomenological purposes, 
it is sufficient to take a simple form for the potential: 
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